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Abstract. Blue compact galaxies (BCGs) are compact 
galaxies that are dominated by intense star formation. 
Comparing the observational properties with the predic- 
tions of stellar population synthesis model, we have an- 
alyzed the nuclear stellar population and emission line 
spectra in a sample of 10 BCGs. The results indicate that 
the continuum flux fractions at 5870A due to old stel- 
lar components and young stellar components are both 
important. The contribution from intermediate age com- 
ponents is different in different galaxies. Our results sug- 
gest that BCGs are old galaxies, in which star formation 
occurs in short intense burst separated by long quiescent 
phases. We have also derived the internal reddening for the 
stellar population by population synthesis method, and 
the internal reddening for the emitting gas clouds by the 
Balmer line ratio. The former is significantly smaller than 
the latter for BCGs. A model of clumpy foreground dust, 
with different covering factors for the gas and stars, can 
explain the difference. Combining the internal reddening 
value and the stellar population, we have decreased the 
effecting from the internal reddening and underlying stel- 
lar absorption, and accurately measured most emission 
lines for each BCGs. Using these emission lines, we have 
attempted to identify the ionizing mechanism of BCGs. 
The ionizing mechanism for these emission line regions of 
BCGs is typical of photoionization by stars, characteristics 
of a low extinction HII regions. 

Key words: galaxies: abundances - galaxies: compact - 
galaxies: star clusters - galaxies: formation - galaxies: stel- 
lar content 



1. Introduction 

Blue compact galaxies (BCGs) were fi rst observed spectro- 
scopically by Sargent & Searle ( 197C ), who clearly estab- 
lished that the properties of these galaxies implied high 
star formation rates at low metallicities (Doublier et al. 
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1997). BCGs have been thought to represent a different 
and extreme environment for star formation compared to 
the Milky Way and many other nearby galaxies. They 
are very important for understanding the star formation 



process and galactic evolution (Kinney et al. 1993, Martin 



1998). BCGs are characterized by their compact morphol- 



ogy and very blue UBV colors (Sage et al. 1992, Hunter 



Thronson 1995). Their optical spectra show strong narrow 



emission lines superposed on a nearly featureless contin- 
uum, similar to the spectrum of HII regions (Izotov et al. 
1997, Ostlin et al. 1999). Radio observations at 21-cm have 



shown BCGs contain large amounts of neutral hydrogen. 
The mean value of Af ni A^ tot for a sample o f 122 BCGs 



is 0.16 (Kriiger et al. |1995i Salzer & Norton |1999D . Sys- 
tematic spectroscopic studies of BCGs have shown that 
about one-third of BCGs have broad W-R bumps, mainly 
at A4650 that are characteristic of late WN stars (Conti 



1991, Izotov et al. 1997) 



Ever since their discovery, the question has arisen 
whether BCGs are truly young systems where star forma- 
tion is occurring for the first time, or whether they are old 
galaxies with current starburst superposed on an old un- 
derlying stellar population (Garnett et al. 1997, Lipovet- 
sky et al. 199£). Because the star formation rate in BCGs 
is very high, the metallicity could have reached the ob- 
served value even within a time T ~ lO^yr (Fanelli et al. 
1988). Hence, one interpretation for the low metallicity, 



high gas content and high star formation rate is that BCGs 
are young objects, and they are being seen at the epoch 
of the formation of the first generation stars (T ~ lO^yr). 
The other interpretation is that they are old objects in 
which star formation occurs in short bursts with long qui- 



escent phases in between (Kriiger et al. 1995, Gondhalekar 



et al. 1995, Ostlin et al. 1999) 



The low metal abundance together with the high star 
formation rates and large gas masses makes BCGs most 
suitable to determine the element abundance (Thuan et 



al. 1995, 1996), the primordial helium abundance Yp (Izo- 
tov et al. 1994) and to study the variations of one chemical 



element relative to another (van Zee et al. 1998). It also 



provides a wealth of diagnostics for the study of intrin- 
sic physical conditions (Izotov & Thuan 1999). The re- 
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suits of these papers are based on direct measurements of 
the emission line intensities, but according to Vaceli et al. 
(1997), the observed emission line intensities are affected 
by the underlying stellar absorption. Since the stellar ab- 
sorption can affect substantially some fundamental emis- 
sion lines used for the derivation of reddening and other 
physical and chemical properties, one of the first and most 
critical steps in the analysis of BCG spectroscopic prop- 
erties is to quantify and remove the contribution of the 
stellar population. 

If we can resolve the stellar population of a BCG, we 
can know its age and star formation regime. We can then 
subtract the stellar absorption line from the emission-line 
spectrum. With the launch of HST and 10-m class tele- 
scope, we are now witnessing a new era that allows to 
analyze in detail nearby objects, such as Galactic HII re- 
gions or 30 Dor in the LMC and resolve old red giants in 
a few distant galaxies (Grebel 1999). Such studies allow us 
to resolve and study individual stars in massive star clus- 
ters. However, as one studies objects at larger distances, 
individual stars (except for some giants) are unresolved 
and hence we are limited to studying their global prop- 
erties (Mas-Hesse & Kunth 1998 ). In this paper we have 
selected 10 BCGs and determined their stellar population 
by applying a population synthesis method based on star 
cluster integrated spectra. In a subsequent paper, we will 
apply an evolution population synthesis method to these 
galaxies and the results of the two papers will be consid- 
ered together. 

The outline of the paper is as follows. In Sect. 2 we 
describe the observations and data reduction. In Sect. 3 
we present measurements of equivalent widths and con- 
tinuum analysis for the BCG spectra. In Sect. 4 we carry 
out the population synthesis and give the results of com- 
putation. In Sect. 5 we subtract the stellar population 
synthesis spectra from the observed ones and study the 
resulting emission line spectra. The results are presented 
in Sect. 6. In Sect. 7, we summarize our conclusions. 
Throughout this paper, we use a Hubble constant of 
Ho = 50 km s^i Mpc~^ 

2. Observations and Data Reduction 

Our sample of 10 blue compact galaxies was selected from 



Kinney et al. 1993; seven of these have Mb > —20 mag 
and are dwarf galaxies (BCDG). Table |l| describes the 
sample properties. The target names are listed in Column 
1. Column 2 to Column 8 respectively give the source coor- 
dinates, morphological type, the radial velocity relative to 
the Local Group, photographic magnitude, absolute mag- 
nitude and Galactic reddening. 

Long-slit spectroscopic observations were carried out 
in March and July of 1997. All the observations were 
made with the 2.16 m telescope at the Xinglong Station 
of Beijing Astronomical Observatory, using a Zeiss uni- 
versal spectrograph with a grating of 300 grooves ■mm~^ 



dispersion. A Tek 1024x1024 CCD was employed, cover- 
ing a spectral range from 3500 to 7500A with a resolution 
of 9.6A (2 pixels). The slit aperture was fixed at 250/iTO, 
corresponding to 2.5" projected on the sky, to match the 
typical seeing value at Xinglong Station, and is set at po- 
sition angle 90° . All the spectra were extracted by adding 
the contributions of 5 pixels around the nucleus. The sky 
background was estimated from a linear interpolation of 
two regions located 30" from the nucleus. The last two 
columns of Table [| list the observation date and exposure 
time. 

The standard reduction to flux units and transforma- 
tion to linear wavelength scale were made using the IRAfQ 
packages. The IRAF packages CCDRED, TWODSPEC 
and ONEDSPEC were used to reduce the long-slit spec- 
tral data. The wavelength was calibrated using a He/Ar 
comparison lamp. 20-odd lines were used to establish the 
wavelength scale, by fitting a first-order cubic spline. The 
accuracy of the wavelength calibration was better than 
2A. On most nights, more than two KPNO standard stars 
were used for relative flux calibration. We followed stan- 
dard procedure in the data reduction: bias subtraction, 
flat fielding, sky subtraction, cosmic ray extinction, CCD 
response curve calibration, wavelength and photometric 
calibration, and extinction correction. The dark counts 
were so low that they were not subtracted. An estimate 
of the radial velocity of the galaxies was made to correct 
to zero red shift. The foreground reddening due to our 
Galaxy was corrected using the value from Burstein & 
Heiles ( 1984 ) . Atmospheric extinction was corrected using 
the mean extinction coefficients for the Xinglong station. 
The final extracted spectra of the BCGs (labeled as OBS) 
are shown in Figure 0. 

3. Measurements 

Our main goal is to study the stellar population of BCGs. 
To do this, we have used the s ynthe sis method described 
in Schmitt, Bica & Pastoriza ( 1991: ). The method mini- 
mizes the differences between the observed and synthetic 
equivalent widths for a set of spectral features. To measure 
the equivalent widths accurately, it is extremely important 
first to have a good fit to the continuum. Hence the analy- 
sis of each of the sample spectra proceeds in two steps (1) 
determining a pseudo-continuum at selected pivot-points, 
and (2) measuring the equivalent widths (EWs) for a set 
of selected spectral lines. 

The continuum and EW measurements followed the 



method outhncd in B ica & AUoin ( |1986| ), Bica (|1988D 
and Bica et al. ( |1994D , Cid Fernandes et al. ( |1998D , and 
subsequently used in several studies of bo th no rmal and 
emission line galaxies (e.g., Jablonka et al. 1990| , Storchi- 
Berg mann et al. 1995 , McQuade et al. 1995 , Bonatto et al. 



1998 ). This method first determines a pseudo-continuum 



IRAF is provided by NOAO 
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Table 1. Parameters of the selected Blue Compact Galaxies and Log of optical observations 



Galaxy 


R.A. 


Decl. 


Morph. 


'^hcljObs 


Bt 


Mb 


E{B-~V)g 


Obs. date 


Exp. 


Name 


(B1950) 


(B1950) 




km/s 


(mag) 


(mag) 




dd/mm/yy 


(sec.) 


IC1586 


00:45:17.0 


22:06:07 




5821 


14.2 


-20.53 


0.02 


25/07/1997 


1500 


NGC2537 


08:09:42.5 


46:08:32 


Sc 


441 


12.4 


-18.13 


0.04 


18/03/1997 


1200 


MRK19 


09:12:53.5 


59:58:53 




4230 


15.6 


-19.07 


0.03 


18/03/1997 


1800 


MRK108 


09:17:26.1 


64:26:57 


10 


1534 


15.4 


-17.57 


0.04 


18/03/1997 


2000 


MRK25 


10:00:22.0 


59:40:50 




2602 


14.8 


-19.96 


0.00 


19/03/1997 


1500 


MRK35 


10:42:16.5 


56:13:23 


Im 


995 


13.2 


-18.50 


0.00 


18/03/1997 


1200 


NGC4194 


12:11:41.7 


54:48:21 


Sm 


2506 


13.0 


-20.71 


0.00 


19/03/1997 


1500 


UGC9560 


14:48:55.1 


35:46:36 


Ir. 


1213 


14.8 


-17.95 


0.00 


18/03/1997 


1500 


UGCA410 


15:35:48.4 


55:25:34 


dE 


665 


15.4 


-16.63 


0.01 


20/03/1997 


2820 


MRK499 


16:47:02.6 


48:47:44 




7710 


14.9 


-21.12 


0.00 


20/03/1997 


2500 



at a few pivot- wavelengths, and then integrates the flux 
difference with respect to this continuum in the deflned 
wavelength windows (Table ^ to determine the EWs. The 
pivot wavelengths used in this work are based on the same 
as those used by the above authors; these were chosen 
to avoid regions of strong emission or absorption features 
(Table |). Four point flux values (3784, 3814, 3866, 3918 
A) were used for the Balmer discontinuity. The use of a 
compatible set of pivot points and wavelength windows is 
important, since it allows a detailed quantitative analysis 
of the stellar population by the synthesis techniques using 
the spectral library of star clusters (Bica & AUoin 1986). 

The determination of the continuum was done inter- 
actively, was taking into account the flux level, noise and 
small uncertainties in wavelength calibration, as well as 
the presence of emission lines. The 5870A point, in par- 
ticular, is sometimes buried underneath the Hel 5876A 
emission line. In such cases, adjacent wavelength regions 
guided the placement of the continuum. The final mea- 
sured fluxes, normalized to the flux at 5870A, are given 
in Table ||. After fitting the continuum points to a con- 
tinuum spectrum, we measured the equivalent widths of 
seven characteristic absorption lines. When a noise spike 
was present in the wavelength window, it was necessary 
to make 'cosmetic' corrections. The results are given in 
Table ^. The first two rows give the names of the spectral 
lines and corresponding wavelength range. The EWs are 
in units of A, and negative sign denotes emission lines. 

4. Stellar Population Synthesis 

To study the history of star formation, age, and ionization 
mechanism of BCGs, we applied the synthesis method of 



Schmitt et al. (1996) for a determination of the stellar pop- 
ulation in their nuclear region. In this method, we start 
with a sample of star clusters, consists of 3 clusters in 
the SMC, 12 in the LMC, 41 Galactic globular clusters 
and 3 rich compact Galactic open cluster, together with 
4 HII regions (Bica & AUoin 1986 ). Then a grid of base 
components (comprising the values of the continuum at se- 



lected points and of the EWs of selected absorption lines) 
is constructed by interpolation and extrapolation in the 
Age — [Z/Zq] plane of those quantities of the clusters of 
the sample. The grid consists of 34 points with ages at 
10^5 X 10"^, 10*, 5 X 10^10^,5 X 10^, and > lO^V, and 
[Z/Zq]^ 0.6, 0.3, 0.0, -0.5, -1.0,-1.5, -2.0, plus one point 
representing the HII region. The method, when applied to 
a given target, consists of adjusting the percentage con- 
tributions of the 35 base components by minimizing the 
difference between the resulting and measured EWs of the 
selected set of absorption lines. When all the resulting 
EWs reproduce those of the galaxy within allowed lim- 
its, we said to have obtained an acceptable solution. We 
then take all the acceptable solutions to form an average 
solution (Schmitt et al. 1996| ). 

The computation can be performed in two ways: 
one way spans the whole Age — [Z/Zq] plane (multi- 
minimization procedure, hereafter MMP), while the other 
is restricted to chemical evolutionary paths through the 
plane (direct combination procedure, hereafter DCP). We 
combine these two methods in our paper. We first use the 
MMP method to single out the main contributing compo- 
nents. And then, based on their resemblance to the whole- 
plane solution and on their reduced chi-square (x^), we se- 
lect the best evolutionary path and use the DCP method 
to give the final result. 



4.1. The results of MMP 

A detailed analysis of the MMP method can be found in 
Schmidt et al. (1989). This method searches the vector 



space of resolutions generated by the entire 35 component 
basis, leading to a representative set of acceptable solu- 
tions to the synthesis problem. We tried various combina- 
tions of the 35 components until a good match between 
the equivalent widths of the synthetic and observed lines 
is obtained. An iterative optimization procedure was used, 
and each iteration alters the percentages of different com- 
ponents. The input parameters are the measured equiva- 
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Table 2. Continuum points Fx relative to Fa(5870A), corrected for the foreground reddening 



Name 




-^3784 


-f3814 




-F3918 


-f4020 


-F4570 


-^4630 


^5313 


-f6630 


IC1586 


1.025 


1.193 


0.996 


1.268 


1.399 


1.311 


1.206 


1.231 


1.045 


0.959 


NGC2537 


0.622 


0.563 


0.739 


0.686 


0.935 


0.891 


0.985 


1.023 


0.907 


0.796 


MRK19 


1.255 


1.189 


1.933 


2.758 


2.201 


1.790 


1.314 


1.249 


0.989 


0.866 


MRK108 


1.108 


0.868 


2.224 


2.207 


1.840 


1.895 


1.041 


1.039 


0.847 


0.886 


MRK25 


1.100 


1.113 


1.215 


1.394 


1.378 


1.338 


1.206 


1.206 


1.089 


0.870 


MRK35 


1.365 


1.078 


1.244 


1.445 


1.559 


1.776 


1.346 


1.338 


1.114 


0.865 


NGC4194 


0.764 


0.776 


0.921 


1.154 


0.957 


1.026 


0.900 


0.922 


0.962 


0.964 


UGC9560 


1.277 


1.538 


1.664 


1.901 


1.843 


1.613 


1.417 


1.295 


0.975 


0.829 


UGCA410 


0.827 


1.422 


1.588 


1.798 


1.540 


1.269 


1.153 


1.114 


1.026 


0.966 


MRK499 


1.012 


1.331 


1.216 


1.442 


1.551 


1.440 


1.252 


1.252 


1.033 


0.935 



Table 3. Measured equivalent widths in A for blue compact galaxies. 



Galaxy 
Name 


KCall 
3908-3952 


4082-4124 


CN 
4150-4214 


4284-4318 


4318-4364 


4846-4884 


Mgl-hMgH 
5156-5196 


IC1586 


5.845 


0.287 


-2.73 


2.122 


-0.89 


-6.90 


1.029 


NGC2537 


18.96 


2.293 


5.692 


3.329 


-14.7 


-33.9 


4.270 


MRK19 


9.878 


5.313 


-0.37 


3.786 


-4.21 


-35.5 


4.138 


MRK108 


18.25 


-8.92 


8.427 


4.192 


-37.1 


-105. 


3.434 


MRK25 


4.318 


3.482 


3.954 


2.682 


-2.39 


-8.96 


2.895 


MRK35 


2.769 


-6.29 


7.278 


1.936 


-15.1 


-22.6 


3.375 


NGC4194 


3.707 


1.261 


0.449 


1.919 


-2.62 


-10.1 


3.226 


UGC9560 


4.736 


-4.67 


4.230 


3.678 


-10.7 


-40.5 


4.303 


UGCA410 


10.61 


1.073 


14.33 


8.688 


-13.0 


-50.6 


3.237 


MRK499 


3.687 


3.092 


-1.07 


0.960 


4.124 


-1.11 


1.161 



lent widths of the selected absorption lines, the continuum 
ratios and a set of trial values of E{B — V) between 0.0 
and 1.0 at steps of 0.02. The result is expressed with the 
flux fractions at 5870 A for each component. The flux frac- 
tions of different components for the 10 BCGs are shown 
in Table ^ (HIIR denotes the HII region component). 

From the results of the MMP method, we find some 
obvious trends for all the BCGs. First, the dominant pop- 
ulation are young (T < 5 x lO^yr) stellar clusters. Second, 
there is a small population of old (T > lO^^yr) and high 
metallicity {[Z/Zq] > 0) globular clusters. Within old 
globular clusters low metallicity components contribute 
more than high metallicity components. Third, for the 
young and intermediate age (T ~ 10^ — 5 x lO^yr) clus- 
ters, components with metallicities below or equal to the 
solar value make a large contribution. It shows that the 
stars in the BCGs have low metallicities, and this is con- 
sistent with our previous knowledge. Lastly, the younger 
components show a large dispersion in the plane, with no 
clear evolutionary paths. This could be due to the fact we 
have only used spectral data in the visible range. Addi- 
tional data in the near ultraviolet and infrared will pro- 
duce better-constrained solutions in the plane. Neverthe- 



less, the presence of starbursts with T w 5 x lO^yr is easily 
recognized from the MMP results. 

4.2. The results of DCP 

To reduce the dispersion in metallicity, all the popula- 
tion synthesis so far made assumed some arbitrarily chosen 
chemical evolution path in the Age — [Z/Zq] plane. The 
improvement in the method of this paper is: we shall pick 
out from the MMP results those components that con- 
tribute most importantly and use them to define paths 
of chemical evolution, thus reducing the degree of arbi- 
trariness. The 3 bright BCGs appear to follow a path 
containing the 11 components along the time sequence 
T = 10^ ~ 5 X lOVr at fixed metallicity [Z/Zq] = -0.5, 
the metallicity sequence [Z/Zq] = —0.5 ~ —2.0 with 
T > lO^Vr, and the HII region. The 7 BCDGs follow 
a path containing the 12 components along the time se- 
quence T = 10^ - 5 X lOVr at fixed metallicity [Z/Zq] = 
0.0, the metallicity sequence [Z/Zq] = 0.0 ~ —2.0 with 
T > 10^°yr, and the HII region. In addition, we also tested 
other possible paths with different maximum metallicities; 
we found that the of the path selected from MMP is 
the smallest. 
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Table 4. Flux fractions (%) at 5870A for the whole plane solution for BCGs. 



5 



HIIR 


E7 


5E7 


E8 


5E8 


E9 


5E9 


> ElO 




HIIR 


E7 


5E7 


E8 


5E8 


E9 


5E9 


> ElO 


a) IC1586 




b) NGC2537 




2.37 


0.45 


0.45 


0.00 


0.00 


0.00 


0.00 


0.6 




1.00 


0.06 


1.36 


2.99 


2.52 


1.23 


0.04 




1.05 


0.45 


0.45 


0.00 


0.00 


0.00 


0.00 


0.3 




0.90 


0.92 


1.16 


3.81 


3.65 


1.06 


0.13 




2.30 


0.45 


0.45 


0.00 


0.00 


0.00 


0.00 


0.0 


0.00 


0.75 


5.79 


5.02 


10.1 


11.5 


6.34 


1.11 


3.36 


6.43 


2.08 


2.04 


14.6 


0.00 


0.00 


0.00 


-0.5 




0.62 


0.16 


0.40 


2.22 


3.04 


2.53 


3.33 












u. / u 


U.Ul 


U.Ul 


-i.U 














^ AP, 
















A 79 


99 


-1.5 














A 


















OU. 1 


-2.0 
















94 








c) MRK19 
















d) MRK108 










0.10 


1.84 


2.33 


8.63 


1.29 


0.00 


0.00 


0.6 




0.00 


0.14 


0.21 


13.5 


1.84 


0.28 


O.GG 




0.10 


1.08 


3.54 


10.3 


1.74 


0.00 


0.00 


0.3 




0.00 


0.00 


0.07 


14.1 


1.90 


1.21 


0.00 


0.00 


0.16 


4.78 


5.97 


26.8 


6.49 


0.13 


0.26 


0.0 


0.00 


0.00 


0.14 


0.28 


31.7 


9.10 


2.13 


0.70 




0.16 


1.78 


2.61 


0.13 


1.82 


0.43 


0.40 


-0.5 




0.00 


0.00 


0.14 


0.00 


1.69 


1.56 


0.78 














1 1 Q 


1 .oo 


-i.U 














9 97 
















1 1 




-1.5 














o.oy 


















A 1 7 


-2.0 
























e) MRK25 
















f) MRK35 










0.16 


0.68 


1.05 


3.74 


2.48 


0.67 


0.02 


0.6 




0.34 


3.33 


3.28 


3.98 


0.27 


0.00 


0.00 




0.17 


1.76 


1.10 


4.68 


2.14 


1.09 


0.34 


0.3 




1.47 


3.10 


2.72 


3.84 


1.99 


0.07 


0.16 


23.7 


0.20 


2.74 


7.11 


10.8 


7.00 


1.71 


0.85 


0.0 


22.2 


4.93 


11.0 


10.6 


12.1 


2.59 


0.97 


0.20 




0.24 


1.66 


2.85 


0.53 


1.31 


1.64 


1.63 


-0.5 




1.23 


2.20 


3.87 


0.00 


0.03 


0.00 


0.00 












z.oo 






-i.U 












u.ou 


n no 

u.uu 


n no 

u.uu 
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O.ZiU 


-1.5 
















u.oy 


















-2.0 
















u.oi 








g) NGC4194 
















h) UGC9560 










4.50 


2.63 


2.64 


0.12 


0.02 


0.03 


0.30 


0.6 




1.90 


1.03 


0.72 


1.82 


1.66 


1.25 


0.27 




4.79 


3.29 


3.40 


0.34 


0.12 


0.11 


0.50 


0.3 




1.23 


1.27 


1.98 


1.73 


1.73 


2.02 


0.27 




5.20 


4.04 


4.04 


0.80 


0.35 
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Table || reports the path solution. The numbers of ac- 
ceptable solutions and corresponding reduced chi-squared 
(X^) are also given. There is a great similarity between 
the path (DCP) and the whole plane (MMP) solution. 
Table 5a, 5g, 5j provide the results for the 3 bright 
BCGs (IC1586, NGC4194 and MRK499). We see that 
the younger components with age 10^ — lO^yr make an 
appreciable contribution, > 35%. The old globular clus- 
ters make even larger contributions. The intermediate-age 
components are small. The other tables in Table || provide 
the results for the 7 BCDGs. The dominant stellar compo- 
nents (> 30%) are in the young age bins (T < 5 x lO^yr). 



The old globular clusters have different values in differ- 
ent galaxies; in some, lower than 10%, and in others as 
much as 35%. The BCDGs differ from the bright BCGs in 
two obvious aspects. First, the intermediate age compo- 
nents in the BCDGs amount to 20%, higher than in the 
bright BCGs. A peak occurs at about 5 x lO^yr, which in- 
dicates that an enhanced star formation event occurred at 
that epoch. Second, the youngest components (T = 10 ''yr) 
in the BCDGs are much lower than in the bright BCGs; 
probably indicating that the star formation rate of BCDGs 
is lower now. The HII region component is a featureless 
continuum, which acts in the synthesis as a dilutor of ab- 
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sorption lines. From Table |[ we note that some BCDGs 
have large contributions from the young components indi- 
cating intense star formation, and small contribution from 
HII regions, which suggests that intense starbursts have 
converted most gas into stars. 

5. Synthesized Spectra 

We now display the results of stellar population synthe- 
sis in a more visible form in Figure g. OBS represents 
the observed spectrum of galaxy (OBS-f 5, the -1-5 indicat- 
ing that it is displaced upwards by 5 units). SYN repre- 
sents the synthetic spectrum resulting from the percentage 
contribution in Table 5. The various stellar components 
are designated as follow. OGC stands for the contribu- 
tion from the old globular cluster (T > lO^^yr); lYC for 
the contribution from the intermediate age star cluster 
(T - 10^ - 5 X lO^r); YBC for the contribution from the 
young blue star clusters 10'' < T < 5 x lO^yr), and HIIR 
stands for the contribution from HII regions. The emis- 
sion line spectrum (OBS-SYN) resulting from subtracting 
SYN from OBS, is shown in the lower part of figure. 

The figure shows that the synthesized spectrum gives 
a good fit to the observed continuum and absorption lines 
for each galaxy. We can conclude that the continuum of 
BCGs comes both from the stars (particularly the young 
and intermediate- age stars) and HII regions. It may be in- 
terpreted that the main energy sources of BCGs are young 
hot O, B stars, which lead to the formation of HII regions 
around them. 

The stellar subtracted spectra (OBS-SYN) can be used 
to study the emission lines. We have measured the main 
strong emission line intensities with Gaussian fits. The re- 
sults, relative to that of Ha are shown in Table At this 
stage our spectra are corrected for the foreground redden- 
ing and the internal reddening from the stellar populations 
(see next section). First, we have used the Hq;/H/3 ratio 
in Table || to derive the internal reddening value associ- 
ated with the line-emitting regions. It will be calculated 
in next section. Second, we have attempted to identify the 
ionizing mechanism in these nuclei, using the emission line 
ratios in the visible region. 

We compared emission-line ratios calculated from 
Table ^ with the diagnostic diagram [NII]6584/Ha — 
[ OII]3 727/[OIII]5007 from Baldwin, Phillips & Terlevich 
( |l98lD . We plot the results in Figure ||. From this fig- 
ure, we find that the BCGs are always located near or 
within the HII region. None of them are located in the 
loci of planetary nebulae, power-law, or shock-heated re- 
gion. This result indicates that the young, massive mass 
stars formed in the nucleus are heating the gas in the nu- 
cleus of BCGs. This result is very similar to those found 
from population synthesis. We also show in this diagram 
the effect of internal reddening. We can see that it, even 
in the case of NGC4194, with the strongest reddening de- 
rived from its line-emitting regions, is not large. 



6. Results 

From the stellar population analysis in Sect. 4 and the 
emission line spectrum in Sect. 5, combined with results 
from other studies, it is now possible to reveal some global 
properties of BCGs. 



6.1. Age and Star Formation Rate 

There are two major competing theories for BCGs. The 
first one claims that BCGs arc truly young systems un- 
dergoing the first star formation episode in the galaxy's 
lifetime. The second model suggests that BCGs are old 
galaxies, which are mainly composed of older stellar pop- 
ulations, with, however, a brief episode of violent star for- 
mation, in order to account for the observed spectroscopic 
features and spectral energy distributions. 

Using the population synthesis method, we know the 
stellar population (different ages and metallicities) per- 
centage contributions at A = 5870A for each of the 10 
BCGs. These results show clearly that, for each, the 
old globular clusters and intermediate age components 
(T < lO^yr) make sizeable contributions to the galactic 
spectrum. The presence of large fractions of old or inter- 
mediate age components indicates that the star formation 
happened already at an early stage, and at a high rate. 
Our results support the second model that BCGs arc old 
galaxies. 

For IC1586, NGC4194 and MRK499, the contributions 
coming from the young and the old stellar components are 
large, but that from the intermediate age component is 
small. It suggests that the rate of star formation during 
the intermediate age period is smaller than in the other 
periods, and the star formation process is not continuous 
in these galaxies. For the BCDGs, the contribution from 
the intermediate age component is important, star forma- 
tion was most vigorous in its intermediate age period, with 
relatively small contributions from the other periods. This 
also implies that star formation in these galaxies is also 
discontinuous. 

The other result of our population synthesis is that 
while the observed properties of the bright BCGs (IC1586, 
NGC4194, MRK499) and the BCDGs are very similar, 
their stellar components and star formation regimes are 
generally different. There are many old and young stellar 
components in bright BCGs, and its recent star formation 
rate is very high. For BCDGs, the old and young stellar 
components are relatively small, but the contribution from 
intermediate age stellar populations is important. 

The stellar population synthesis suggests that BCGs 
are old galaxies, in which the process of star formation 
is intermittent. Star formation has been violent in one 
of its evolution periods. These results are also supported 



by ot her observations (Papaderos et al. 1996, Sung et al 



1998 , Aloisi et al. 1999). It illustrates that the present 



method is more than a simple population synthesis since 
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Fig. 1. Observed spectra for our 10 BCGs are shown together with the best-match population synthesis models (cf. 
Table H). The observed galaxy spectrum (OBS) has been corrected for the foreground reddening E(B-V)g' and shifted 
by 5 units. In the synthetic spectrum (SYN), the young population (YBC) dominates the light, but the intermediate 
age component (lYC) is also important. The emission line spectrum appears in the (OBS-SYN) difference, in the lower 
part of the figure. OBS and SYN are normalized to F> = 10 at 5870A. 
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Table 5. Flux fractions (%) at 5870A for the best path solution of BCGs. 
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it provides a direct estimate of the chemical evolution of 
the galaxy. 



6.2. Internal Reddening 

When we investigate the internal energy source, physi- 
cal condition and internal structure of galaxies, we must 
take into account the effect of internal reddening (Piza- 
gno & Rix 1998). The effect of dust extinction on the 



emerging radiation is one o f the least unde rstood physical 
phenomena (Calzetti 1997 , Ho et al. 1997). To study the 



internal reddening properties of BCGs, we quantify the 



discrepancy between the dust extinction measured from 
the emission line ratios and the optical continuum. 

In the method of population synthesis used in this pa- 
per, the internal reddening is taken as an adjustable pa- 
rameter, so that an estimate for the internal reddening is 
made at the same time as the stellar composition. We try 
various values of the internal reddening, make the appro- 
priate correction in the continuum spectrum, then use this 
corrected continuum spectrum in the synthesis to find the 
best solution. This is an empirical way of determining the 
galactic reddening, its advantage is that it is assumption- 
free. The values of galactic internal reddening are listed 
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Table 6. Emission line intensities with respect to Ha (normalized to 100) and the internal reddening value of BCGs. 
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Fig.l continue. 



in Table ||, E{B - V)mmp is the result from the MMP 
method, E{B — V)dcp is the result from the DCP method. 
We find that the values are small {E{B—V) < 0.35), which 
is consistent with the BCGs being metal-poor and dust- 
poor. We find the reddening also clearly depends on the 
shape of the spectrum. The flattest spectrum (NGC4194) 
goes with the largest color excess. The steeper spectrum, 
the less the extinction. 

The Balmer line ratio Ha/H/3 allows us to character- 
ize the dust extinction in the regions where the nebular 
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Fig. 2. Logarithmic plot of the line intensity ra- 
tios [NII]A6584/Ha versus [Off] A3727/[OIII]A5007. Open 
symbols correspond to data corrected for the foreground 
reddening. Black symbols represent data points with an 
additional reddening correction intrinsic to the Hne emit- 
ting regions. 

lines are produced. We measured the internal reddening 
value E{B — y)HQ/H/3 of the 10 BCGs using the observed 
emission lines Ha and H/?. The difference in the calcula- 
tion of E{B — l^)Ha/H/3 between the previous work and 
ours is that we can correct the underlying stellar absorp- 
tion EWafcs from the results of stellar population synthesis 
without making any hypotheses. The result of E(B — V) 
is listed in the last 3 columns of Table ^. 

From this table, we can find that the internal reddening 
of the stellar continuum in BCGs is generally lower than 
that of ionized gas. A model of foreground dust clumps, 
with different covering factors for gas and stars, is a pos- 
sible explanation for the difference. The covering factor 
by dusty clumps is greater for the gas region that gener- 
ates the emission lines than for the stars that produce the 
continuum. That the continuum emission of stars is less 
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obscured than are the emission hnes of ionized gas, has 
been pointed out for other kinds of emission Une galaxies 



(Calzetti et al. 1994) 



7. Summary 

We have observed the optical spectra for the nuclear re- 
gions of 10 blue compact galaxies. We have studied their 
stellar populations by matching the spectra of these ob- 
jects to a library of integrated spectra of star clusters. Our 
conclusions can be summarized as follows. 

The quantitative analysis indicates that the nucleus of 
the 3 bright BCGs is dominated by young components and 
the star-forming process is still ongoing. The maximum 
metallicity of the stellar population is [Z/Zq] = —0.5. 
The nucleus of the dwarf BCGs (BCDGs, which have 
Mb > —20), on the other hand, is dominated by the in- 
termediate age component. The metallicity has, at most, 
reached up to the solar value. The young component is 
not so important as in the bright BCGs, but it is still not 
negligible. 

For all BCGs, the old population in the range 
\ZIZq\ < 0.0 is important, the very metal-rich compo- 
nent provides quite a small contribution in these galaxies. 
The stellar populations of BCGs suggest that they are old 
galaxies with intermittent star formation history. 

A good match can be achieved between the synthesized 
and observed spectrum of BCGs. It suggests that the stel- 
lar radiation is an important energy source for BCGs. 

The emission-line spectra from the gaseous compo- 
nents in these objects were isolated and analyzed. Using 
these ^tellar siibtra.cted spectra,, we have calculated the 
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